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The  study  of  phase  changes  in  overlayers  is  an  important  tool  in 
understanding  the  interactions  that  adatoms  undergo. Very  little  work 
has  been  done  on  chemisorption  systems,  and  only  H/W(100)^^^  and  0/W(110)^^^ 

have  been  studied  in  any  detail.  This  Letter  reports  a new  phase  transition 

t 

for  0/W(110)  at  low  coverages.  We  observe  what  we  believe  to  be  two- 

I 

I dimensional  dissolution  of  the  chemisorbed  islands  previously  identified  in 

■ (3  41 

this  system'  ’ As  a result,  we  are  able  to  determine  separately  attractive 
i and  repulsive  adatom-adatom  interaction  energies  leading  to  the  p(2  x 1) 

I 

structure. 

(31 

The  experiment'  ' consists  of  measuring  the  superlattice  beam  intensities 

and  angular  profiles  in  low-energy  electron  diffraction  (LEED)  from  the  W(llO) 

surface  covered  with  varying  amounts  of  oxygen.  The  most  important  feature  of 

this  system  is  that  the  ordered  overlayer  forms  by  an  island  growth  mechanism. 

This  is  established  by  the  observation  that  superlattice  LEED  beams  form  already 

(4l 

at  very  low  coverage'  . It  implies  a net  attractive  adatom-adatom  interaction. 

However,  the  p(2  x 1)  structure  that  is  formed,  shown  in  Fig.  1 and  consisting 

of  doubly  spaced  close-packed  rows  parallel  to  <111>  directions,  requires  a 

nearest-neighbor  repulsion.  Streaking  of  spots  is  not  observed  in  the 

(3  4a  5l 

diffraction  pattern  during  formation  of  the  overlayer'  ’ ’ . This  implies 

that  island  formation  is  preferred  over  the  formation  of  long  rows,  and 
, therefore  that  the  attractive  interactions  along  <111>  directions  are  about 

equal  in  magnitude^^’®^  and  considerably  smaller  than  the  nearest-neighbor 
repulsion  that  prevents  close  packing. 

I 


Oxygen  coverage  determinations  were  made  by  establishing  saturation 
coverage  for  the  p(2  x 1)  structure  on  the  basis  of  maximum  superlattice  beam 
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intensity  and  using  sticking  coefficient  measurements^^^  to  relate  exposure 
time  to  relative  coverage.  The  temperature  dependence  of  superlattice  beams 
at  a given  coverage  was  then  measured  for  various  diffraction  goemetries. 

In  each  case  the  data  plotted  on  a Inl  vs  T graph  consist  of  a region  of 
constant  slope  proportional  to  the  Debye-Waller  factor  for  thermal  vibrations 
in  the  direction  of  the  diffraction  vector,  and  a sudden  falloff  away  from 
this  line.  Figure  2 shows  examples  of  this  falloff  (with  the  Debye-Waller 
contribution  subtracted),  for  a low  coverage  of  oxygen  and  for  the  p(2  x 1) 
saturation  coverage.  At  the  low  coverage,  the  superlattice  beam  intensity 
begins  to  decay  at  a much  lower  temperature.  Figure  3 shows  the  transition 
temperature  as  a function  of  coverage,  where  is  chosen  as  the  inflection 
point  in  the  curves  of  Fig.  2.  Clearly  there  are  two  transitions,  one 
operative  at  low  coverage,  the  other  for  saturation.  In  the  intermediate- 
coverage  region,  both  transitions  are  seen,  the  lower-coverage  one  becoming 
weaker  as  the  coverage  increases,  but  remaining  at  the  same  T^. 

We  have  interpreted  the  low-coverage  transition  as  the  dissolution  of 
ordered  regions,  with  atoms  leaving  the  ordered  regions  and  distributing 
themselves  randomly  on  sites  in  the  "sea"  or  uncovered  part  of  the  W substrate. 
To  do  this,  only  the  attractive  interaction  leading  to  formation  of  ordered 
regions  needs  to  be  overcome,  and  in  that  sense  the  transition  is  similar  to 
the  sublimation  of  a solid  or  dissolution  of  a solute  into  solvent  and 
fundamentally  different  from  the  saturation-coverage  transition,  where 
disordering  can  occur  only  by  atoms  moving  into  the  more  strongly  repulsive 
nearest-neighbor  sites.  The  latter  thus  corresponds  to  a true  order-disorder 
transition,  with  a higher  transition  temperature,  as  shown  in  Fig.  3. 


-4- 


Since  the  dissolution  and  order-disorder  transition  temperatures  are 

respectively  a measure  of  the  attractive  interaction  leading  to  island 

formation  and  the  repulsive  interaction  preventing  close  packing,  it  now 

becomes  possible  to  determine  these  interactions  separately  by  first  fitting 

the  low-coverage  transition  and  then  using  the  resulting  attractive  interaction 

to  fit  the  saturation  - coverage  transition  and  obtain  the  repulsive  interaction. 

Because  of  the  overlayer  symmetry,  one  assumption  is  still  necessary,  i.e. 

that  the  attractive  interactions  - and  - (see  Fig.  1)  are  equal. 

To  model  the  low-coverage  transition,  we  consider  the  dissolution  of 

islands  consisting  of  an  r-sublattice  on  which  all  sites  are  occupied  and 

f8l 

an  s-sublattice  on  which  all  sites  are  empty'  . As  the  temperature  increases, 

the  total  number  of  atoms  in  the  island,  and  thus  participating  in  the  diffraction, 

changes.  Hence  this  is  a problem  in  the  grand  canonical  ensemble,  and  the 

(91 

lattice  gas  model  developed  by  Lee  and  Yang‘  ' gives  a satisfactory  descrTption 

of  the  system.  The  present  phase  transition  then  is  the  two-dimensional  analog 

(91 

of  the  condensation  problem,  to  which  the  Lee  and  Yang  model'  ' has  been 

successfully  appl ied^^^’^ ^ \ Defining  b = 1 if  an  r-sublattice  site  in  the 

island  is  occupied,  and  b = 0 if  an  r-sublattice  site  is  empty,  or  for 

"l"2 

any  site  in  the  sea,  correspondence  can  be  made  to  the  ferromagnetic  model 
for  the  spin  system  in  zero  external  magnetic  field^^^^  by  setting 

n ' n ~ 1 = 1 if  an  r-sublattice  is  occupied  /-.x 

"l”2  ^l"2  = -1  if  an  r-sublattice  is  empty  ' ' 
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where  are  integers  to  indicate  lattice  positions, 

operator  for  site  (n^n^)  and  -J  is  the  usual  nearest-neighbor  attractive  spin- 
spin  interaction  energy  in  the  Ising  model.  Then 


<b 


^1^2  ^1^2 


b . .>  = 1/4  [<a, 


a,  ,>+<a  >+<a,  i>+l]  (3) 

n^jn^  n^n^  ■' 


^1^2  ^1^2 


is  the  pair  correlation  function  with  <a  the  usual  spin-spin 

”2  ^1 ^2 


correlation  function  and  <a  _ > the  average  magnetization  in  the  spin  system. 

"r2 

As  is  well  known,  the  angular  distribution  of  intensity  in  any  diffraction 

1)2) 

experiment  is  a measure  of  the  pair  correlation  function'  , the  two  being 
related  by  a Fourier  transform.  In  our  case  the  LEED  superlattice  reflections 
are  a result  of  the  new  periodicity  introduced  by  the  ordered  regions  of  0 atoms. 
Hence  the  angular  distribution  of  intensity  in  these  beams  is  a measure  of 
the  overlayer  pair  correlation  function.  Evaluating  this  at  the  diffraction 
maximum  gives  for  the  maximum  LEED  intensity  scattered  by  the  overlayer 


c c 


■max'"’  ' . ^ J..  "■’n.n,  <">  '’nin'  1^)’ 


^1^2  ^1^2 


r’2 


(4) 


where  the  f's  are  structure  factors  that  are  functions  of  the  adsorbed  layer 

island  size  and  take  into  account  the  effects  of  boundaries  on  the  diffracted 

intensity,  and  is  the  coherence  area  of  the  instrument. 

We  have  calculated  the  correlation  functions  of  Eq.  (3)  for  different 

Ca  = e2b‘  evaluated  I„,g^(T)  through  Eq.  (4),  and  then  compared  with 

experimental  r,„(T)  curves.  A typical  fit  is  shown  in  Fig.  2a.  The  detailed 
rriox 

shape  of  the  model  intensity  decay  with  temperature  depends  both  on  the 
correlation  functions  and  the  f. 's  of  Eq.  (4),  but  e.  does  not.  For  arbitrary 

1 d 


I • 
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finite  lattice  size,  the  correlation  functions  are  very  complicated  and  no 

(131 

analytic  forms  are  known'  ' . We  have  estimated  these  functions  using  the 
(141 

Onsager  solutions'  ' for  an  infinite  two-dimensional  lattice.  Details  of 

(151 

this  are  presented  elsewhere'  . The  f.'s  in  Eq.  (4)  are  the  most  important 
quantity  in  determining  the  shape  of  An  experimentally  determined^^^ 

^ iTIaX 

0 

-35A  dia  mean  island  size  was  used  to  fit  the  curves.  Fits  to  the  shape  of 

Imax(T)  with  assumed  50%  greater  or  lesser  island  sizes  were  much  poorer. 

(14l 

The  interaction  energy  e,  was  determined  from  the  Onsager  relation' 

a 

2 El  2 E. 

sinh  -r-  sinh  7 i-  = 1 . (5) 

"^B't  ^B't 

where  here  Ei=E2=l/4  e^.  A fit  to  the  experimental  transition  temperature  of 

460°K  gives  £3=^26  eV/atom.  As  expected  from  Eq.  (5)  the  determined 

value  of  is  a quite  sensitive  function  of  T^. 

This  value  for  the  attractive  interaction  can  now  be  used  to  analyze 

the  second  phase  transition,  for  the  saturation  p(2xl)  coverage,  T^.  = 718°K. 

No  further  assumptions  are  necessary  to  extract  the  repulsive  interaction 

that  prevents  the  formation  of  a close  packed  structure.  However,  for 

this  phase  transition,  a different  model  is  required  than  for  the  dissolution 

transition,  because  the  number  of  atoms  participating  in  the  diffraction 

is  now  constant.  The  model  we  use  to  describe  it  is  a generalization  of 
(21 

a model'  ' considering  isotropic  n.-n.  interactions  only.  In  our  case 
anisotropic  n.-n.  interactions  -e^  and  +C|^  are  required.  This  is  equivalent 
to  a ferromagnetic  system  in  zero  magnetic  field  with  attractive  n.-n. 
interactions  -J^  and  if  the  following  correspondence  is  made: 
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-^a  • "S  i 

With  this  generalization,  we  use  the  model  of  Ref.  2 to  calculate 

again  the  correlation  functions  and  I , (T),  the  maximum  intensity  of  the  | 

nicix  I 

^ 1 

superlattice  reflection  as  a function  of  T,  in  similar  fashion  as  we  did  Ij 

earlier  for  the  dissolution  phase  transition.  We  fit  the  experimental  H 

j 

^max^^^  curves  for  the  saturation  p(2xl)  coverage  by  calculating  correlation  j 

functions  for  different  e.  using  e,  from  our  earlier  results.  As  before  the 

H A i 

transition  temperature  is  defined  through  Eq.  (5).  This  gives  = 0.15  i 

■ to) 

eV/atom.  This  value  lies  between  that  obtained  by  a more  approximate  model'  ' 
and  that  obtained  by  an  assumed  interaction  strength  in  a Monte-Carlo 
calculation^^’^^^ 

It  seems  evident  that  dissolution  transitions  should  be  observable  for 
other  island-forming  adsorbates.  Since  the  fit  of  the  above  models  to  T^ 
does  not  depend  critically  on  instrument  response  or  island  size,  it  should 
be  straightforward  to  extract  adatom  interaction  energies  for  such  adsorbates. 


I. 
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FIGURE  CAPTIONS 


t Figure  1.  Structure  model  of  p(2  x 1)  island  on  a bcc  (110)  substrate. 

I 

X's  indicate  sites  occupied  by  0 atoms,  labeled  the  r-sublattice 
j in  the  text.  Filled  circles  indicate  unoccupied  sites.  Within 

an  island,  these  are  termed  the  s-sublattice.  ^^>^2b  '^®P'^®sent 
interactions  between  sites  as  indicated  by  the  unit  vectors. 

Figure  2.  Temperature  dependence  of  the  (1/2  1/2)  superlattice  reflection 
at  high  and  low  coverage  and  comparison  to  model  calculation. 

a)  low  coverage  6 = 0.25; 

b)  saturation  p(2  x 1)  coverage  0 = 0.53.  Note  the  diffe>'-''nt 
transition  temperatures  and  shapes  of  the  curves. 

Figure  3.  Experimental  transition  temperature  vs.  coverage.  0 = 0.5 

corresponds  to  saturation  coverage  for  the  p(2  x 1)  structure. 
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